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ABSTRACT 

We perform numerical simulations of impulsively generated Alfven waves in an isolated photospheric 
flux tube, and explore the propagation of these waves along such magnetic structure that extends 
from the photosphere, where these waves are triggered, to the solar corona, and analyze resulting 
magnetic shells. Our model of the solar atmosphere is constructed by adopting the temperature 
distribution based on the semi-empirical model and specifying the curved magnetic field lines that 
constitute the magnetic flux tube which is rooted in the solar photosphere. The evolution of the solar 
atmosphere is described by 3D, ideal magnetohydrodynamic equations that are numerically solved by 
the FLASH code. Our numerical simulations reveal, based on the physical properties of the multi¬ 
shell magnetic twisters and the amount of energy and momentum associated with them, that these 
multi-shell magnetic twisters may be responsible for the observed heating of the lower solar corona 
and for the formation of solar wind. Moreover, it is likely that the existence of these twisters can be 
verified by high-resolution observations. 

Subject headings: Sun: atmosphere. Sun: corona, magnetohydrodynamics (MHD), magnetic fields 


1. INTRODUCTION 

The existence of hot solar corona (> 10^ K) is well- 
established by observations of solar ultraviolet and X- 
ray emissions originating in magnetically confined coro¬ 
nal plasma, and the acceleration of fast solar wind that 
commence from magnetically open coronal regions (Pe¬ 
ter and Dwivedi 2014). The identification of physical 
processes responsible for the coronal heating and solar 
wind acceleration has been a very challenging and still 
remains unsolved problem (Priest et al. 1998, Cranmer 
2002, Peter and Dwivedi 2014). Recent high-resolution 
solar observations seem to clearly imply that the energy 
and mass transport required to heat the corona and ac¬ 
celerate the solar wind must come from the chromosphere 
through highly localized wave and dynamical plasma pro¬ 
cesses at small spatio-temporal scales, such as spicules, 
network jets, swirls, and twisting motions of magnetized 
plasma (De Pontieu et al. 2004, Shibata et al. 2007, De 
Pontieu et al. 2014, McIntosh et al. 2011). 

Solar chromospheric swirls have been detected by re¬ 
cent observations (Wedemeyer-Bohm et al. 2012). They 
are categorized by their various morphological shapes 
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(Tian et al. 2014). These swirls are suggested to be 
driven by solar convective motions (Wedemeyer- 
Bohm et al. 2012, Steiner et al. 2010, Shelyag 

et al. 2011). Magnetically-confined plasma in 
such swirls couples to the solar upper atmospheric lay¬ 
ers, transferring swirls rotation into the low corona. 
At higher altitudes, these swirls may appear like large 
cyclones/tornadoes due to expansion of magnetic flux 
tubes that support them (Shelyag et al. 2011, Su et 
al. 2012, Zhang and Liu 2011). The relationship be¬ 
tween the chromospheric swirls and coronal tornadoes 
is not yet understood (Wedemeyer et al. 2013, Shelyag 
et al. 2013, Wedemeyer and Steiner 2014). It is also 
debated whether the swirls are triggered by solar con¬ 
vective motions, or they are associated with torsional 
Alfven waves at small spatial scales (Shelyag et al. 2013, 
Wedemeyer and Steiner 2014). The observations, how¬ 
ever, suggest that these localized motions of magnetized 
plasma might transfer up to 1.4 x 10^ Wm“^ en¬ 
ergy flux to the base of the solar corona for its localized 
heating and wind acceleration (Wedemeyer-Bohm et al. 
2012). 

The main objective of this paper is to describe novel so¬ 
lar phenomena called here multi-shell magnetic twisters, 
and their driving mechanism and related morphology. It 
should be noted that recently observed chromospheric 
swirls can be an integrated appearance of multi-shell 
magnetic twisters. Our finding sheds new light on the 
morphological evolution of such twisters, their exact ex- 
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citation mechanism and their role in the localized heating 
of the solar corona, and mass transport to nascent solar 
wind. The novel result is our theoretical prediction of 
the multi-shell magnetic twisters, whose existence can be 
verified by observations with upcoming high-resolution 
instruments. 

The paper is organized as follows. Our model of so¬ 
lar magnetic arcades and description of our numerical 
method are introduced in Sects. 2 and 3, respectively. 
Results of our numerical simulations of the Alfven wave 
propagation in magnetic arcades are presented and dis¬ 
cussed in Sect. 4. Conclusions are given in Sect. 5. 

2. NUMERICAL MODEL OF MULTI-SHELLS 
2.1. MHD equations 

We consider a gravitationally stratified and magneti¬ 
cally confined plasma in a structure that resembles a flux 
tube, which is described by the following set of MHD 
equations: 

If+ V-(eV) = 0,(l) 

^ (V • V) V =-Vp + 1(V X B) X B + pg ,{2) 
at fi 

— = Vx(VxB),(3) 
V-B = 0,(4) 
^+V-Vp = -7PV-V,(5) 

p=^pT,( 6 ) 

m 

where g is mass density, p gas pressure, V, B and 
g = ( 0 ,—^, 0 ) represent the plasma velocity, the mag¬ 
netic field and gravitational acceleration, respectively. In 
addition, T is temperature, m particle mass, that was 
specified by mean molecular weight value of 1.24 (Oskar 
Steiner, private communication), is the Boltzmann’s 
constant, 7 = 1.4 the adiabatic index, and p the mag¬ 
netic permeability of plasma. The value of g is 274 m 
s“^. 

2 . 2 . A model of the statie solar atmosphere 

We consider a model of the static {d/dt = 0) solar 
atmosphere in which all plasma quantities vary with x, 
and z. We assume that the solar atmosphere is in static 
equilibrium (Vg = 0) with the Lorentz force balanced 
by the gravity force and the gas pressure gradient, which 
means that 

1(V X Be) X Be + Qeg - Vpe = 0 , (7) 

P 

where the subscript 'e' corresponds to the equilibrium 
configuration. We consider an axisymmetric flux tube of 
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Fig. 1.— Magnetic field lines at t = 0 s. 
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Fig. 2.— Numerical blocks used in the numerical simulations 
and the initial pulse, V 0 {x, y, z = 0), drawn in the vertical {x — 
y) plane for z = 0. Note that a part of the whole simulation 
region is displayed only and max{\Ve\) = 2.56 km . 


which equilibrium is described by Murawski et al. (2015) 
which is based on the analytical theory developed by 
Solov’ev (2010). This flux tube is initially non-twisted, 
and its magnetic field is represented by the azimuthal 
component of the magnetic flux function {A 0) as 

Be = Vx(M^), 


(8) 
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where ^ is a unit vector along the azimuthal direction. 
In this case, we have 


B 


er 


dA p p _ 1 d{rA) 

, ^e9 — U , ney — -^- 

oy r or 


(9) 


where r = is the radius. For a flux tube, we 

can specify A{r, y) as 

A{r, y) = Bo exp {-k^y"^) ^ > (10) 

where Byo is the magnitude of the external magnetic 
field which is chosen vertical, and ky are inverse 
length scales along the radial and vertical directions, re¬ 
spectively. The magnetic field lines, which follow from 
Eqs. (|9]) and (pd)]) are displayed in Fig. [T] for k^ = ky = 4 
Mm“^. We choose and hold fixed the magnitude of the 
reference magnetic field Bq in such way that the magnetic 
field within the flux tube, at r = ^ = 0 Mm, is about 
285 Gauss, and Byo ^11.4 Gauss. For these values, the 
resulting magnetic field lines are predominantly vertical 
around the tube axis, r = 0 Mm, while further out they 
are bent and Bq decays with distance. Such magnetic 
field lines correspond to an isolated axi-symmetric mag¬ 
netic flux tube. It must be also noted that the magnetic 
field at the top of the simulation region is essentially uni¬ 
form, with its value of ByQ. 

Having specified the magnetic field, the equilibrium 
mass density and the gas pressure are evaluated from 
Eq. ([7]) with the adaptation of the hydrostatic gas pres¬ 
sure ph which is given by 


Ph{y) =Po exp 


r dy 

Jyr Hy') 


( 11 ) 


where y^ = 10 Mm is the reference level, po = 0-01 Pa is 
the gas pressure evaluated at y = y^, noting that 


Hy) 


ksThiy) 

mg 


( 12 ) 


is the pressure scale-height, and T\^{y) a hydrostatic tem¬ 
perature profile. We adopt T\^{y) for the solar atmo¬ 
sphere that is specified by the model developed by Avrett 
& Loeser (2008). This temperature profile is smoothly 
extended into the corona. It should be further noted 
that in our model the solar photosphere occupies the re¬ 
gion 0 < p < 0.5 Mm, the chromosphere is sandwiched 
between p = 0.5 Mm and the transition region located 
at y 2.1 Mm, and above this height the atmospheric 
layers represent the solar corona. 


3. NUMERICAL SIMULATIONS OF MHD EQUATIONS 

In order to solve Eqs ([I])-® numerically, we use the 
ELASH code (Eryxell et al. 2000; Lee & Deane 2009; 
Lee 2013), in which a third-order unsplit Godunov-type 


solver with various slope limiters and Riemann solvers 
as well as Adaptive Mesh Refinement (AMR) (MacNeice 
et al. 1999) are implemented. The minmod slope lim¬ 
iter and the Roe Riemann solver (e.g., Toth 2000) are 
used. We set the simulation box as ( — 1.5 Mm, 1.5 Mm) x 
(0 Mm, 18 Mm) x (—1.5 Mm, 1.5 Mm) and impose fixed in 
time boundary conditions for all plasma quantities in all 
directions. In our present work, we use a static, non- 
uniform grid with a minimum (maximum) level of re¬ 
finement set to 2 (5). Note that small size blocks of 
numerical grid occupy the level up to p = 4 Mm (Eig. [2j), 
and every numerical block consists of 8 x 8 x 8 identical 
numerical cells. This results in the smallest spatial reso¬ 
lution of 23.4 km below the level p = 4 Mm and allows 
to well resolve spatial structures in the low solar corona. 


3.1. Initial perturbations 


The atmospheric magnetic field is continuously dis¬ 
turbed by large-scale dynamical perturbations that are 
able to transfer kinetic energy {e.g. buffeting due to 
the granular motion in the photosphere, or flare-driven 
blast waves in the upper atmosphere, reconnection- 
driven shocks, etc.). We assume that the instigator of 
changes in the magnetic field has such a nature. How¬ 
ever, its exact property is not specified, and initially (at 
t = 0 s) we perturb the above described equilibrium im¬ 
pulsively by a Gaussian pulse in the azimuthal compo¬ 
nent of velocity, V^, viz.. 


Ve = Ay— exp 
w 


r'^ + {y- yof 




(13) 


where = {x — x^Y 4-z^ is the squared radius. Ay is the 
amplitude of the pulse, (xo,po,0) its position, and w its 
width. We set Ay = 6 km s“^, w = 150 km, and po = 500 
km, and hold them fixed, while we allow xq to attain 
one of these two values: (a) xq = 0 km; (b) xq = 100 
km. This value of Ay results in the effective maximum 
velocity of about 2.56 km s“^ (Eig. [SJ color map), and 
the value of po shows that the system is perturbed right 
at the top of the photosphere. 

Note that in our 3D model, the torsional Alfven waves 
decouple linearly from magnetoacoustic waves. They can 
be described solely by Ve{x^yA) (see Murawski et al. 
2015). As a result, the initial pulse triggers Alfven waves. 


4. RESULTS OF NUMERICAL SIMULATIONS 

We simulate small amplitude and impulsively excited 
Alfven waves and investigate their propagation along 
the magnetic field lines which are associated with the 
flux tube (Eig. [T]). We consider the following two cases: 
(a) centrally-launched initial pulse; (b) off-centrally- 
launched initial pulse, and describe the corresponding 
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Fig. 3. — Spatial profiles of magnetic field lines at t = 200 s 
for the case of centrally-launched initial pulse. 

results in the following parts of the paper. 

4.1. Centrally launched pulse 

The case of a centrally-launched pulse corresponds to 
the location of the initial pulse at the flux tube axis, 
which is realized by setting xq = 0 Mm in Eq. m- A 
3D view of the magnetic field lines of the flux tube is 
shown at t = 200 s in Fig. [H It is clear that as a result 
of torsional Alfven propagation, magnetic field lines are 
twisted and this twist is localized essentially below the 
transition region. 

Figure |4] presents velocity streamlines which exhibit 
the twisters’ motion. At t = 290 s we discern a centrally 
located shell in which plasma counter-rotates in compar¬ 
ison to the rotation of the external shell (the top panel). 
At the later moment of time, this central shell is more 
developed (the bottom panel). It is to be further noted 
that these concentric magnetic shells have their dimen¬ 
sion of the order of 120 km. 

Figure [5] shows the vertical cross-sections of 
Vz{x,y,z = 0) = Ve{x^y^z = 0) exhibiting the ori¬ 
gin of various magnetic interfaces (multi-shells) at an 
upper chromospheric height of ^ = 1.9 Mm, which is 
located about 200 km below the transition region. It 
is to be noted that the magnetic field perturbations 
are confined to atmospheric layers located below the 
solar transition region (Murawski and Musielak 2010; 
Murawski et al. 2015), while the velocity perturba¬ 
tions easily penetrate into the corona and transfer 
the available mechanical energy (Fig. [5]). The vertical 
cut at the upper chromospheric height in the velocity 
clearly shows the two oppositely oriented shells, where 
plasma is counter-streaming to each other and forming 
two twisting plasma shells. As a result, a distinct set 
of magnetic interfaces is formed, where Alfven wave 
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Fig. 4.— Spatial profiles of streamlines at t = 290 s (top) 
and t = 420 s (bottom) for the case of centrally-launched 
initial pulse. 

velocity perturbations are also generated. 

Figure [6] illustrates spatial profiles of the vertical com¬ 
ponent of velocity, Iy(x, y^z = 0, t), drawn in the vertical 
plane for 2 ; = 0. As a result of the initial pulse in V^, it 
follows from the Bernoulli’s law that the under-pressure 
region is generated shortly after t = 0 s and the ambient 
plasma is sucked into the region where the initial pulse 
was launched from. This under-pressure sucks the ambi¬ 
ent plasma, resulting in a converging into the launching 
region flow. Such flow is well seen at t = 50 s (the top- 
left panel). Note the down-flowing plasma in the patch 
located at (x = 0,^ ^ 0.75) Mm, while the up-flowing 
gas is located at (x = 0,^ ^ 0.3) Mm. Later on, at 
t = 100 s these patches move upward and new regions of 
up-flowing plasma are generated at {x ~ ±0.2,^ 1.4) 

Mm (top-right). Some of the up-flowing plasma is at- 
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Fig. 5. — Spatial profiles of transverse velocity Vz{x,y,z = 
0, t) at t = 290 s (top) and t = 420 s (bottom) for the case of 
centrally-launched pulse, xo = 0 km. 

tracted by the gravity force, resulting in a down-flowing 
blob of the plasma velocity of ~ —6 km s“^, which 
is seen at t = 290 s above the transition region in the 
form of the violet map (the bottom-left panel). This 
free-falling plasma clashes with the first jet which orig¬ 
inates right below the transition region (the red patch 
in the bottom-left panel). As a result of this clash a 


shock is generated. The second jet is seen at t = 420 s 
(bottom-right). 

The vertical components of velocity, 4^, are displayed 
in the horizontal plane just below the transition region, 
at ^ = 1.9 Mm, in Fig.[3 The development of concentric 
shells is well seen at t = 290 s (top); the inner shell is 
represented by the violet ring with down-flowing plasma 
with its speed Vy ^ 3 km s“^, while in the outer shell, 
marked by the red ring, and the plasma flows upward 
with its speed Vy ^ 1 km s“^. In the center the plasma 
begins to flow upwards. At t = 420 s this central flow 
develops the structure with 4 discernible arms (bottom). 
Note that the vertical flow is out of phase to Vz at each 
instance in a particular shell at a given chromospheric 
height (compare with Fig. [5j). This is a unique corre¬ 
lation of periodic changes in the vertical flows (i.e., up- 
and down-flows) negatively coupled with the Alfven wave 
velocity perturbations in a chromospheric shell. 

The multi-shells can also be traced on the horizon¬ 
tal profiles of mass density, g{x,y = 1.9, z,t), which are 
displayed in Fig. [51 The initial profile (not shown) is 
altered at t = 290 s with 2 well-developed dense concen¬ 
tric shells (top). At a later moment of time, t = 420 
s, there are more shells seen (bottom). These shells re¬ 
sult from up-flowing and down-flowing plasma and from 
the magnetoacoustic-gravity waves, which propagate in 
the horizontal direction. These magnetoacoustic waves 
are surface-gravity-like waves and they are driven by the 
oscillations in Vy . These waves propagate along the tran¬ 
sition region. The propagating upward plasma brings a 
dense gas from lower atmospheric regions, while down¬ 
falling flow leads to plasma evacuation, which is well seen 
in the rings. Within the center of the flux-tube, the os¬ 
cillations in Vy{x,y = 1.9,2;) and g{x,y = 1.9,2;) are in 
anti-phase. Compare Figs. Hand [51 

Figure [9] (top) displays time-signatures of the transver¬ 
sal (Vz) and vertical (Vy) components of velocity, repre¬ 
sented respectively by solid and dashed lines, and col¬ 
lected at the detection point which is located at (x = 0.1 
Mm, ^ = 1.9 Mm, 2 ; = 0.1 Mm). These time-signatures 
reveal perturbations in which Vz and Vy are approx¬ 
imately in anti-phase. Hence, we infer that up-flows 
(which correspond to 1^ > 0) lead to clock-wise oscil¬ 
lations in Ve (represented in the vertical plane by K), 
while down-flows are associated with anti-clock-wise os¬ 
cillations in Vq. The Alfven wave energy flux can be 
approximated by the following relation (Vigeesh et al. 
2012 ): 

Ef « Q.Vg^CA , (14) 

where g^ is the unperturbed mass density, V6» is the az- 
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the case of centrally-launched pulse, xo = 0 Mm. 
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Fig. 7. — Spatial profiles of Vy{x,y = 1.9, 2:, t) at t = 290 s 
(top) and t = 420 s (bottom) for the case of the centrally- 
launched initial pulse, xo = 0 km. 

imuthal component of the Alfven speed and ca is the 
Alfven speed which is expressed as 


c\{x,y,z) 


Bl{x,y,z) 
liQe{x,y,z) ' 


(15) 


Note that the vertical flows associated with the magnetic 
twister carry the mass-flux of its maximum of 6 x 10“^ 
kg m“^ s“^ to the inner coronal heights, and thereby 
they also contribute to the mass supply in the region of 
nascent solar wind formation (Figure [H middle). The net 


Fig. 8. — Spatial profiles of mass density g(x,y = 1.9, z,t) 
at t = 290 s (top) and t = 420 s (bottom) for the case of 
centrally-launched initial pulse. 

energy flux carried by these twisters is about 0.2 x 10^ W 
m“^ (Fig. [HI bottom), which is almost sufficient for the 
required coronal heating of the quiet-Sun and the solar 
wind energy losses. 

4.2. Off-centrally-launched pulse 

We now consider the case of off-centrally-launched 
pulse that corresponds to the location of the initial pulse 
out off the flux tube axis, which is realized by setting 
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Fig. 9. — Time-signatures of Vz (top, solid line), Vy (top, 
dashed line), vertical component of the mass density flux gVy 
(middle), and approximated Alfven waves energy flux (bot¬ 
tom), evaluated at the detection point (x = 0.1 Mm, y = 1.9 
Mm, z = 0.1 Mm) for the case of xo = 0 km. 
xq = 0.1 Mm (Eq. [13]). The spatial profiles of stream¬ 
lines are drawn in Fig. [TOl At t = 150 s, we observe 
a development of the inner and outer eddies, as well as 
side vortices, seen at the bottom-left and top-right cor¬ 
ners (top). At a later moment of time, t = 300 s the 
flow pattern is more complex with more local eddies gen¬ 
erated as a result of energy cascade into smaller scales 
(bottom). By comparing this with Fig. 2] we conclude 
that the off-central pulse results in a more complex flow 
than the central pulse. This conclusion is further sup¬ 
ported by the results of Fig. [11] that illustrates vertical 
profiles of Vz{x,y,z = 0) and Fig. [12] which shows the 
horizontal profiles of g{x,y = 1.9, z). The symmetry of 
the centrally-launched pulse (Fig. [5|) is broken for the 
case of the off-centrally-launched pulse (Fig. [11]), with 
complex structures developed in the latter case. 

5. SUMMARY AND CONCLUSIONS 



Fig. 10. — Spatial profiles of streamlines at t = 150 s (top) 
and t = 300 s (bottom) for the case of xo = 100 km. 


In this paper, we simulated impulsively generated, ei¬ 
ther centrally or off-centrally, Alfven waves in a gravita¬ 
tionally stratified and magnetically confined solar mag¬ 
netic flux tube with the temperature profile of Avrett and 
Loeser (2008). A novel result is that an outer shell with 
Alfven wave velocity perturbation is generated within 
the flux tube at small spatial scales. The plasma associ¬ 
ated with this shell rotates following the magnetic field 
lines at this interface with the velocity streamlines form¬ 
ing the swirls. The vertical flows are associated with this 
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Fig. 12.— Spatial profiles of mass density g{x,y = 1.9, 2 :, t) 
at t = 150 s (top) and t = 300 s (bottom) for xo = 0.1 Mm. 


Fig. 11.— Spatial profiles of transverse velocity Vz(x,y,z = 
0, t) at t = 150 s (top) and t = 300 s (bottom) for the case of 
xq = 100 km. 

dynamics and they feedback energy to generate one more 
inner shell with Alfven wave velocity perturbation within 
the tube. This process repeats in time resulting in a 
third inner shell. On a time-scale of few hundred seconds, 
the multi-shell magnetic twister is fully developed in the 
tube. When we view it from top of the tube, the col¬ 


lective appearance shows up as the perturbations in the 
magnetic field lines with the streamline rotational mo¬ 
tion constituting the multi-shell magnetic twister. Each 
distinct appearance of the multiple plasma shells grows 
in time as every shell is demarcated by the transition 
in density, which results in the emission / ^ N^. The 
scenario of multi-shells development is more complex in 
the case of the off-centrally launched Alfven pulse, with 
a larger number of asymmetric shells generated in the 
system. 
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Fig. 13. — The line-of-sight view of spatially averaged mass 
density-map in the units of kg m“^ (a) at the upper 

chromospheric height below the transition region is compared 
to the observational features of recently discovered ring-type 
chromospheric swirl (Wedemeyer et al. 2013) (b). The com¬ 
parison shows that recently observed chromospheric swirls 
may be an integrated appearance of our multi-shell magnetic 
twisters; note that the spatial resolution of the presented ob¬ 
servations is 540 km per pixel. 


The main new result is that the multi-shell magnetic 
twister and associated multiple torsional Alfven waves 
are generated within a magnetic flux tube embedded in 
the solar atmosphere. The net energy flux carried by 
these twisters is about 0.2 x 10^ W m“^. This is al¬ 
most sufficient for the required coronal heating of the 
quiet-Sun and the solar wind energy losses. The verti¬ 
cal flows associated with the magnetic twister carry the 
mass-flux of its maximum of 6 x 10“^ kg m“^ s“^ to the 
inner coronal heights, and thereby they also contribute 
to the mass supply in the region of nascent solar wind 
formation. The Poynting energy flux carried out by the 
magnetic twisters is high compared to the energy trans¬ 
ported by chromospheric swirls. Therefore, we suggest 
that the magnetic twisters must be considered as an im¬ 
portant candidate for the localized coronal heating and 
the solar wind acceleration. 

The shell solutions seem to be the natural consequence 
of centering the initial pulse on the axis of the tube. This 
ensures that the perturbations have the same symmetry 
as the tube, and rings/shells are the obvious result. For 
more realistic perturbations, where the symmetry is not 
exact, multiple vortices with different orientations are 
likely to be the result (see e.g. the ’realistically’ driven 
vortices in Moll et al. 2012, their Figs. 9 and 11). For 
the implications of this on the observed heating see van 
Ballegooijen et al. (2011). 

Based on the current observations, the typical mor¬ 
phologies of the swirls are: ring and spiral (Wedemeyer 
et al. 2013). If we fade the resolution of density map 
as shown in Fig. 8 by 2-3 times, then the multi-shell 
chromospheric swirls appear as a ring type swirl as re¬ 


cently observed in the solar chromosphere (Wedemeyer 
Bohm & Rouppe van der Voort, 2009; Wedemeyer-Bohm 
et al. 2012, Wedemeyer et al. 2013). The ring-type 
swirl and its apparent motion are intimately connected 
with the localized magnetic flux tube, filled-in plasma, 
and their collective dynamics under the influence of the 
driver. The multi-shell magnetic twister will appear like 
ring-type swirl in a straight flux tube when observed with 
imager with less spatial resolution. This comparison is 
shown in Fig. [131 where an observed ring-type chromo¬ 
spheric swirl matches in morphology over almost similar 
spatial-scale of modeled faded (density averaged) multi¬ 
shell magnetic twister. However, our model twister will 
appear as multi-shell magnetic twister once the spatial 
resolution is scaled down to less than 70 km per pixel. 

Therefore, our newly proposed model of multi-shell 
magnetic twister provides the exact underlying physics 
of their drivers, morphology, and forecast on observa¬ 
tions. The comparison shows that recently observed 
chromospheric swirls may be an integrated appearance of 
our multi-shell magnetic twisters. Moreover, the present 
novel model solves the paradigm of the association of 
Alfven waves with such twisters. It is shown that they 
are present at multi-shells and collectively carry large 
amount of energy to heat the corona locally and acceler¬ 
ate the solar wind. 

Based on the physical properties of the multi-shell 
magnetic twisters and the amount of energy and mass 
carried by them, it is suggested that these multi-shell 
twisters are responsible for the observed heating of the 
solar inner corona and for the formation of nascent so¬ 
lar wind. Moreover, it is likely that the existence of the 
multi-shell magnetic twisters can be verified by observa¬ 
tions at higher resolutions. 
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